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ABSTRACT 



We present a comprehensive theoretical study, using a semi-analytical model within the standard LCDM framework, 
of the photometric properties of the progenitors of present-day early-type galaxies in the redshift range < z < 1. 
We explore progenitors of all morphologies and study their characteristics as a function of the luminosity and local 
environment of the early-type remnant at z = 0. In agreement with previous studies, we find that, while larger 
early-types are generally assembled later, their luminosity-weighted stellar ages are typically older. In dense cluster-like 
environments, ~ 70 percent of early-type systems are 'in place' by z = 1 and evolve without interactions thereafter, while 
in the field the corresponding value is ~ 30 percent. Averaging across all environments at z ~ 1, less than 50 percent of 
the stellar mass which ends up in early- types today is actually in early- type progenitors at this redshift, in agreement 
with recent observational work. The corresponding value is ~ 65 percent in clusters, due to faster morphological 
evolution in such dense environments. We develop probabilistic prescriptions which provide a means of including spiral 
(i.e. non early-type) progenitors at intermediate and high redshifts, based on their luminosity and optical colours. 
For example, we find that, at intermediate redshifts (z ~ 0.5), large (My < —21.5), red (B — V > 0.7) spirals have 
~ 75 — 95 percent chance of being an early-type progenitor, while the corresponding probability for large blue spirals 
(Mb < —21.5, B — V < 0.7) is ~ 50 — 75 percent. The prescriptions developed here can be used to address, from 
the perspective of the standard model, the issue of 'progenitor bias', whereby the exclusion of late- type progenitors in 
observational studies can lead to inaccurate conclusions regarding the evolution of the early-type population over cosmic 
time. Finally, we explore the correspondence between the true 'progenitor set' of the present-day early-type population 
- defined as the set of all galaxies that are progenitors of present-day early-types regardless of their morphologies - and 
the frequently used 'red-sequence', defined as the set of galaxies within the part of the colour-magnitude space which 
is dominated by early-type objects. We find that, while more massive members (My < —21) of the 'red sequence' 
trace the progenitor set reasonably well, the relationship breaks down at fainter luminosities (My > —21). Thus, while 
the results of recent observational studies which exploit the red sequence are valid (since they are largely restricted to 
massive galaxies), more care should be taken when deeper observations (which will probe fainter luminosities) become 
available in the future. 



Key words, galaxies: elliptical and lenticular, cD 
parameters 



galaxies: evolution - galaxies: formation - galaxies: fundamental 



1. Introduction 

As 'end points' of galaxy merger sequences, early-type 
galaxies carry important signatures of mass assembly and 
star formation in the Universe. Deducing their star forma- 
tion histories (SFHs) therefore contains the key to under- 
standing not only the evolution of these galaxies but the 
evolutionary patterns of galaxies as a whole. Our view of 
early-type galaxy formation has developed over the years, 
awa y from the classical 'monolithic col lapse' hypothesis 
(e.g. lLarsoi]|1975HChiosi &: Carraroll2002|) and towards the 
hierarchical assembly of these objects through mergers and 
accretion of smaller galaxies over time, in the framework 
of the cu rrently popular LCDM paradigm of galaxy forma- 
tion (e.g. Kauffmann et al.lll993l,ll996al iBaugh et al.l 1996 ; 



Cole et all 120001: lHatton et alJl2003t iDe Lucia et "all 120061 



Bower et al.l 20061 and references therein 



iKauffmann fc Charlotl 119981: ISomerville fc Primackl [L999 : 



skaviraj@astro.ox.ac.uk; s.kaviraj@imperial.ac.uk 



A significant body of observational work in the past 
has traced the assembly histories of early-type galax- 
ies by stud ying only early-type populations at high 
redshift (e.g iGladders et al.l [l99l (Stanford et alJ 119981: 
iBarrientos fc Lillvl l2003t iFerreras et al.ll2005D . However, a 
fundamental feature of early-type formation in the stan- 
dard LCDM model is that stellar mass that eventually 
ends up in present-day early-type galaxies is not entirely 
contained in early-type systems at high redshift. Although, 
early-types at any redshift are almost guaranteed progen- 
itors of their counterparts at present-day, looking only at 
early-types at high redshift introduces a 'progenitor bias', 
which becomes increasingly more severe at larger look-back 
times, as the fraction of early-type galaxies becomes pro- 
gressively smaller and late-type sy stems begin to domi- 
nate the progenitor population (e.g. iFranx fc van Dokkuml 
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19961: Ivan Dokkum et al.1l2000Hvan Dokkum fe Fra^l200lt 
Kavirai et al.ll2005l) . 

In the current era of large scale surveys, e.g. SDSS 
(Adelman-McCarthy et al. 2006), COMBO-17 (Wolf et al. 
2004), MUSYC (Gawiser et al. 2006), GEMS (Rix et al. 
2004), unprecedented amounts of data spanning a large 
range in rcdshift (typically < z < 1) and environment 
are becoming available, allowing us to study statistically 
significant numbers of galaxies at various stages of evolu- 
tion. A quantitative study of early-type progenitors within 
the standard model is therefore desirable to (a) understand 
the difference in the way early-types are assembled as a 
function of their luminosity and environment and (b) to 
gauge the role of non-early-type progenitors in early-type 
evolution, especially for studies that focus at high redshift. 

A central theme of this work is to use a model which 
accurately reproduces the photometric properties of early- 
type galaxies and their observed evolution over a large 
range in redshift to study, in detail, the properties of the 
progenitors of early-type galaxies at present day. The pro- 
cess of mapping the progenitor population and its evolution 
with redshift also provides a realistic picture of progenitor 
bias with in the framework of the stand ard model. 

While van Dokkum fc Franxl (|200lD have developed and 
studied the issue of progenitor bias, their study employed 
phenomenological SFHs, with the simple (and perhaps un- 
realistic) assumption that morphological transformations 
occur abruptly ~ 1.5 Gyrs after the cessation of star forma- 
tion in a particula r galaxy. Their work was an extension o f 
previous ideas fe.g lBower et al.lll998HShiova fc Be"kk"Il 998') 
that took into account the potential for elliptical galaxies to 
have complex star formation histories, but did not explore 
the effect of morphological transformations on the evolution 
of the early- type population at high redshift. 



While some recent 


observational studies 


(e.s. iHolden et al. 20051: 


van de Ven et all 120031: 


van Dokkum & van der Marel 


20071) have used the 


results of van Dokkum & Franx (2001) to take the effects 



of progenitor bias into account, other studies (e.g. Blakeslee 
et al. 2003) have used the entire galaxy population, without 
reference to morphology, in an attempt to reduce the bias 
in their conclusions that would occur if only elliptical 
galaxies were used in the analysis. 

Our study refin e s an d extends the results of 
Ivan Dokkum fc Franxl (|200ll ) by studying early-type pro- 
genitors within a realistic and well-calibrated semi- 
analytical framework, in which mass assembly and mor- 
phological transformations can be followed more accu- 
rately in the context of the LCDM paradigm. While the 
Ivan Dokkum fc Franxl (|2001f ) study does attempt to allevi- 
ate the effects of progenitor bias, its simplicity, especially 
in the prescription used for morphological transformations, 
makes a more realistic treatment of this issue, in the frame- 
work of the standard model, very desirable. Using the en- 
tire galaxy population without reference to morphology is 
also not ideal because that implies that all galaxies at high 
redshift are potential progenitors of present-day early-type 
galaxies. While this approximation may be reasonably ro- 
bust for the most massive early- type galaxies in the densest 
regions of the Universe, our study shows that this is gener- 
ally not applicable. 

In this paper we study the evolution of the 'progenitor 
set' of early- type galaxies, defined as the set of all galax- 
ies at a given redshift that are, regardless of morphology, 



the progenitors of an early- type at z = 0. Note that 'early- 
types' are defined as galaxies that are elliptical or lenticular. 
See Section 2 for a description of how galaxy morphologies 
are defined in the model. We explore the evolution of the 
progenitor set with redshift, as a function of the luminosity 
and environment of the early-type remnant which is left at 
present-day. We pay particular attention to spiral progen- 
itors in the model, since these might have been excluded 
from some studies of early-type evolution in the past, by 
virtue of their morphology, even though they form an im- 
portant part of the progenitor set. By comparing the prop- 
erties (optical colours and luminosities) of spiral progeni- 
tors to the general spiral population, we provide a means 
of correcting for progenitor bias, which is consistent with 
the properties of the standard model, by including specific 
parts of the spiral population at high redshift into the study 
of early-type evolution. 

The plan of this paper is as follows. In Section 2, we de- 
scribe the salient features of the model used in this study. 
In Section 3, we quantify the morphologies of the galaxies 
that make up the progenitor set. We map the properties of 
elliptical, SO and spiral progenitors as a function of redshift 
and explore differences between progenitors as a function 
of the mass and environment of the elliptical remnant at 
present day. In Section 4 we focus exclusively on spiral pro- 
genitors and compare their photometric properties to the 
general spiral population, to derive probabilistic prescrip- 
tions for including spiral progenitors in early-type studies. 
Section 5 traces the contribution of galaxies in dense re- 
gions at high rcdshift to cluster early-types at present-day. 
Finally, in Section 6 we explore the correspondence between 
the true progenitor set of present-day early-type galaxies 
and the 'red-sequence', which is sometimes used as a proxy 
for the progenitor set in observational studies (e.g. Bell et 
al. 2004, Fabcr et al. 2007). 

Note that throughout this study we provide rest-frame 
magnitudes for all model galaxies. Unless otherwise noted, 
the filters used are in the standard Johnson system. 



2. The model 

The semi-analytical model used in this study is GALICS, 
which combines large-scale cosmological N-body simula- 
tions with analytical recipe s for the evolution of baryons 
within dark matter haloes (jHatton et al.ll2003h . GALICS 
makes predictions for the overall statistical properties of 
galaxy populations, with an emphasis on their panchro- 
matic spectral energy distributions across a wide wave- 
length range (UV to infrared/submillimetre). In this sec- 
tion we describe the salient features of the model, includ- 
ing the implementation of the cosmological simulation and 
details of the semi-analytics that determine the growth 
of bulges. We also highlight the successes of GALICS in 
terms of reproducing the photometric properties of early- 
type galaxies in the £/U-optical spectral ranges across a 
wide redshift range. 

Before we begin our analysis we briefly describe the def- 
inition of galaxy morphology in the model used in this 
study. Galaxy morphology in the model is determined 
by the ratio of the B-band luminosities of the disk and 
bulge components which correlates well with Hubble type 
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(jSimien fc de Vaucouleursl Il986h . A morphology index is 
defined as 



exp 



Li 



(1) 



such that a pure disk has I = 1 and a pure bulge has / = 
0. Following lUaueh et all (|l996l ). ellipticals have I < 0.219, 
SOs have 0.219 < / < 0.507 and spirals have / > 0.507. 

2.1. Dark matter simulation 

The dark matter 'backbone' used in GALICS has been gen- 
erated by simulating a LCDM model in a cube with a co- 
moving size of 100/i _1 Mpc with parameters f2 m = 0.33, 
= 0.667 and er 8 = 0.88. The amplitude of the power 
spectrum is computed using the present-day abundance of 
rich clusters (|Eke et al. 199Q) and initial conditions a re ex- 
tracted using the GRAFIC code (jBertschingerlll995l ). The 
simulation contains 256 3 dark matter particles. The mass 
of each dark matter particle is ~ 8 x 10 9 M Q and the spatial 
resolution achieved is ~ 29.3 kpc. The minimum mass of 
galaxies (due to the minimum number of particles used in 
the friend-of-friends group finder and the baryon fraction 
used) is - 2 x 10 10 M Q . 

2.2. The growth of bulges in the semi-analytical model 
2.2.1. Mergers vs disk instabilities 

GALICS includes two processes which lead to the formation 
and growth of spheroidal components in galaxies: mergers 
(both major and minor) and gravitational disk instabilities 
Q. While mergers directly reveal the hierarchical nature of 
galaxy formation, disk instabilities can result in the for- 
mation of bulges in a 'monolithic' manner (i.e. in objects 
classified as isolated at our resolution limit). The relative 
contribution of each of these two processes in the build-up 
of bright, local galaxy bulges thus yields a natural estimate 
of how hierarchical galaxy formation truly is. 

For each of our galaxies we measure the stellar and gas 
masses transfered from the disk to the bulge during merg- 
ers {M merg ) and/or disk instabilities (M inst ). In Figure [TJ 
we show the relative contributions of mergers (both major 
and minor) and disk instabilities to the mass build-up of 
the spheroidal components of local galaxies brighter than 
M(B) = —18.9. This is shown as the distribution of the 

ratio M merg /M instl computed for each model galaxy. A 2.2.3. Post-merger morphology 



mergers to other processes in the formation of spheroids. 
The fact that model elliptical galaxies develop their bulges 
mainly through mergers rather than through disk instabil- 
ities is an output of our model - we do not simply assume 
that elliptical galaxies only form during mergers but in- 
stead allow the different physical processes we believe are 
relevant to the formation process to compete. In that sense, 
one can consider the results we obtain as a demonstration 
that the formation of ellipticals is indeed hierarchical. 

Second, elliptical and lenticular galaxy formation ap- 
pears to be somewhat different, since GALICS predicts that 
most lenticular galaxies develop their bulges through disk 
instabilities. Changing the modelling of disk instabilities 
would affect the plots in a similar way for ellipticals and 
lenticulars, leaving the difference in their respective be- 
haviours unchanged. The fact that the model lenticulars 
acquire their morphology through disk instabilities is a re- 
sult which is common in numer ical simulations of galactic 
dynamics (|Combes et al. 2000) and gives us some confi- 
dence that our simple physical modelling broadly captures 
the mechanism. 

Finally, the lower panels of Figure Q] indicate that spi- 
rals form their bulge components mainly through disk in- 
stabilities. The global picture for morphogenesis that we 
obtain from GALICS is satisfactory as it compares well to 
more detailed numerical simulations. Moreover, we show 
that whilst an 'isolated object' approximation can be jus- 
tified to describe spiral galaxies (and some lenticulars), it 
cannot be used to properly model the formation of the vast 
majority of ellipticals in the framework of CDM structure 
formation. 

2.2.2. A note on the treatment of disk instabilities 

The disk instability model employed is similar to that of van 
den Bosch (1998). Rotational equilibrium is enforced at the 
disk half mass radius - in other words, the circular velocity 
Vtot arising from the presence of total mass of gas, stars 
and dark matter enclosed in the half mass radius sphere 
is calculated. It is then compared to the circular velocity 
V c of the disk alone, and if this latter is smaller than a 
critical value (obtained by mass weighting the value of 0.7 
for a pure gas disk and 0.52 for a pure stellar disk), the 
minimum amount of mass of gas and stars required to re- 
establish stability is transferred to the disk. 



galaxy with a value of 1 for this ratio has acquired half 
its bulge mass during mergers and the other half through 
disk instabilities. Note that we only consider galaxies with 
a bulge component. Pure disk galaxies (i.e. those that have 
never been unstable or undergone a merger) represent only 
~ 1.4 percent of the bright population. 

Several important conclusions can be drawn from this 
figure. First, it is clear that the majority of elliptical galax- 
ies develop their bulges during mergers rather than disk 
instabilities. This result may seem natural, as GALICS is 
an implementation of the hierarchical galaxy formation sce- 
nario. However, the main point here is that, for the first 
time, we do indeed compare the relative contribution of 



In the literature mergers are typically modelled by tak- 
ing the ratio of the progenitors, adding the stars of 
the lighter galaxy to the disk of the heavier one if 
the mass ratio is less than fbulge ~ 0.3 or destroy 
the disk and form a bulge if the ratio is higher (e.g. 
Wal ker et aT1ll996t ICole et al II200QI: ISomerville fc Primackl 
1999; Kauff mann et all I1999H Since a galactic disk can 
be disrupted by an encounter even when the interloper is 
less massive than the disk itself, this simple prescription 
reproduces this behaviour. However, the sharp cut-off at 
fbulge ~ 0.3 between totally disrupting or not disrupting 
the morphology of the galaxy seems somewhat unrealistic. 



1 A disk is considered to be stable if V c < 0.7Vt o t, where V c is 
the circular velocity of the disk and Vtot is the rotational velocity 
of the disk-bulge-halo system (e.g. Ivan den Boscr][l99Sl ). 



In keeping with the literature, we define a merger as 'major' 
when the mass ratio between merging objects ranges between 
one and a third, and as 'minor' when it is smaller. 
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Fig. 1. The ratio between stellar and gas masses transferred 
from the disk to the bulge during mergers (M merg ) and 
through disk instabilities (Mi nst ). A galaxy with a value 
of 1 for this ratio has acquired half its bulge mass during 
mergers and the other half through disk instabilities. This 
plot is restricted to model galaxies brighter than M(B) = 
— 18.9. Note that the (artificial) peaks at 4 and -4 are for 
ratios M merg /Mi nst = oo and respectively. The vertical 
dashed lines represent median values. 



The model implemented in GALICS is constructed in 
terms of a smooth function (X) , which models the fraction 
of disk material remaining in the disk after the merger as a 
function of the mass ratio of the two progenitors. The previ- 
ous studies (described above) effectively use a step function, 
in the sense that the material kept in the disk is either or 
1 depending on whether the mass ratio is greater or lesser 
than 0.3 respectively. In GALICS this smooth function is 
defined as: 



X(R) 



1 



R- 1 



(2) 



where R represents the mass ratio of the heavier to the 
lighter progenitor and x is the critical value of that ratio i.e. 
the value that X — 0.5. Since step functions with fbuige ~ 
0.3 have been found to give good results in the past, the 
fiduc ial value adop t ed in GALICS is x = 1/0.3. See Figure 
3 in lHatton et al.l (|2003| ) for a visualisation of X{R) for 
this fiducial v alue o f x- We refer readers to Section 5 in 
lHatton et al.l (|2003F) for details of prescriptions related to 
merging employed by GALICS. 

Galaxies are modelled with three components: the, disk, 
the bulge and the burst. The burst has the same geometry 
(and therefore star formation law) as the bulge but its scale 



radius is 10% of that of the bulge. During a merger, a frac- 
tion of X of both gas and stars originally in the disk remain 
in the disk while the rest are transferred to the burst. Stars 
formed/existing in the burst are rapidly transferred to the 
bulge. As the central burst forms stars (typically over short 
timescales than the 'quiescent' mode due to its smaller scale 
radius), the new stellar component gets transferred to the 
bulge. As it subsequently runs out of the gas deposited by 
the merger, the burst eventually disappears. We refer read- 
ers to section 5 of Hatton et al. (2003) for further details. 
Note that, for higher values of R (i.e. mergers with a high 
mass ratio) almost all disk material remains in the disk, 
while equal mass mergers result in rapid bulge formation, 
mimicking the prescriptions employed in the literature. 

2.3. A discussion of free parameters in the model 

Inevitably, any study that exploits a semi-analytical model 
such as the one presented here, relies on simplified, albeit 
well-calibrated recipes to describe galaxy evolution. These 
semi-analytic recipes are driven by free parameters which 
determine the evolution of the baryonic Universe. While the 
parameters that drive the baryonic evolution in the fiducial 
GALI CS model are discussed in Section 7 of lHatton et al.l 
(2003), we briefly revisit them here and discuss their po- 
tential impact on the analysis that follows in the paper. In 
Table 1 we list the free parameters in the model and in- 
dicate how their values are constrained in this study. The 
first section of the table indicates parameters that have the 
greatest impact on the evolution of the galaxy population 
and the bottom section indicates parameters that play a 
minor or negligible role. 

While the model is driven by several free parame- 
ters, those that have the greatest impact on the spectro- 
photometric properties of the galaxy population are (a) the 
star formation efficiency (b) the efficiency of mass loading 
for supernova feedback and (c) the merging law that deter- 
mines morphological transformations during galaxy inter- 
actions. We explore each of these in turn and briefly discuss 
the sensitivity of the model to these parameters. 

Star formation efficiency: The star formation effi- 
ciency (SFE) determines the rate at which cold gas is 
converted into stars. For example, decreasing the SFE 
results in less cold gas being converted into stars, moves 
the peak of the star formation rate closer to present 
day and affects both the shape and normalisation of 
the z=0 galaxy luminosity function. In our fiducial 
model, the inverse of the SF E (0) is set equal to 50 (i.e. 
the SFE is -2%), following iKennicutd (| 19981 see their 
Eqn. 7), who combined Ha, HI, CO and far-infrared 
measurements in spiral and infrared-selected starburst 
galaxies to provide a parametrisation of the global 
star formation rate in the local Universe. Figure 6 
in their study indicates that this parametrisation is 
accurate over a wide gas density range, from gas-poor 
spiral disks to the cores of the most luminous starburst 
galaxies, implying that the value for this parameter 
that is used in the model is reasonably well-constrained. 
Furthermore, the present-day B and K-band galaxy 
luminosity functions predicted by GALICS are very 
consistent with those observed by the 2dF survey 
(Cross et al. 2001), indicating that the effect of the 
adopted value (/? = 50) is consistent with the observed 
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distribution of galaxy luminosities in the real Universe. 

Mass loading for supernova feedback: For a given 
stellar mass formed, a certain fraction is contained 
in massive, fast-evolving stars that become super- 
novae (SN), injecting kinetic energy into the interstellar 
medium. The efficiency of this SN dri ven 'wind' depends 
on both the porosity of the fSM (see Silk 2001) and the 
'mass loading factor', which describes the entrainmen t 
of interstellar gas by the wind (see e.g. [Si^[200lll2003l) . 
For example, increasing the mass loading efficiency will 
produce more feedback, heat more cold gas and reduce 
the amount of fuel available for star formation. This 
will, in turn, affect th e galaxy lum i nosity function at 
present-day. Following iMartin et al.l (|2002f ). we assume 
that the mass loading factor is ~10, which results in 
the outflow rate of a starburst to be of the order of the 
star formation rate (see Section 4.2 in Hatton et al. 
(2003) for more details). Thus the inverse of the mass 
loading factor (e) is set to 0.1. As mentioned before, 
the reproduction of the galaxy luminosity functions in 
the B and K bands - which relies in part on e - are well 
reproduced by the fiducial GALICS model. 

Merger law: A third parameter that has a significant 
impact on galaxy evolution is the merger law that 
determines post-merger morphology. This law deter- 
mines the morphological mix of the Universe at z — 
and therefore has a direct impact on the early-type 
progenitor set in the model. The prescription used to 
determine post-merger morphology has already been 
descri bed in detail above in Section 2.2.3. IHatton et aLI 
(2003) indicates that the predicted morphological mix 
in GALICS is (within counting errors) consistent with 
that seen in the real Universe. Note, how ever, that the 
compa risons presented in Section 8.6 of IHatton et al.l 
(2003) use reasonably small surveys (e.g. the Stromlo- 
APM redshift survey; Loveday et al. 1996) compared to 
the scale of modern surveys such as the SDSS. Future 
papers will present more robust comparisons between 
the morphological mix predicted by GALICS and that 
observed in the SDSS e.g. through comparison of the 
GALICS predictions to visually-inspected morphologies 
of the e ntire SDSS DR6 me asured by the 'Galaxy Zoo' 
project (jLintott et al.ll2008Tl . 

Minor parameters: We complete our description of the 
free parameters by briefly discussing minor parameters 
that do not produce a measurable impact on the anal- 
ysis. The normalization for the satellite-satellit e merg - 
ing law (tp) is constrained using Mak ino fc Hull ([1997). 
Satellite-satellite merging is much rarer than dynamical 
friction merging (i.e. merging with centrals in a halo) 
and thus this parameter plays a negligible role in the 
analysis. 

The recycling efficiency (?) describes the fraction of gas 
originally expelled from the halo by feedback which is 
re-accrete d over time as the DM halo gr ows in mass. We 
note that ISomerville fc Primackl (|1999h set this param- 
eter to (i.e. an y gas expel l ed is lost from the system 
altogether) while lCole et al.l (|2000D set it to 1. Hatton et 
al. (2003) use a fiducial value of 30%, to represent an in- 
termediate behaviour compared to the two extremes de- 
scribed above, which is plausibly closer to the truth. We 



note, however, that the effect of this parameter is very 
weak because at the mass resolution of the model, very 
little gas is expelled from DM halos. The re-accreted gas 
is always significantly less than the amount of 'pristine' 
gas accreted or the amount of gas that is ejected out of 
the galaxy and into the halo in the first place. 
In GALICS star formation driven by disk instabilities 
and mergers takes place in a central 'burst' region, 
which is modelled as having the same morphology as 
the bulge but with a fraction of its radius (see Section 4 
in Hatton et al. 2003 for a more details). The 'burst to 
bulge' radius (k) - set to a fiducial value of 0.1 - dictates 
the size of the burst region in model galaxies (i.e. an or- 
der of magnitude smaller than the radius of the bulge). 
Reducing k would make the burst regions smaller, short- 
ening the star formation timescales (and vice- versa). 
However, typical star formation timescales (a few tens of 
Myrs at most) in the burst regions are negligible com- 
pared to the timescales over which galaxy properties 
are being 'viewed' in the model (several Gyrs). Hence, 
changing this parameter does not affect the observed 
colours of the galaxy population predicted by the model 
and thus leaves our conclusions unchanged. 
The dust recipe used in the model is calibrated using 
the Milky Way, the Large Magellanic Cloud and Small 
Magellanic Cloud and a few local spirals for which the 
extinction curv e, gas content and metallicity have been 
measured (see iGuiderdoni fc Rocca-Volmerangd I1987L 
for details). While a key assumption is that the dust 
properties of galaxies are invariant with redshift, tests 
of this assumption requires data at high redshift which 
are not yet available, leaving us very little room to 
further calibrate our du st recipe. Finally , the fiducial 
GALICS model uses the iKennicuttl (jl983h Initial Mass 
Function (IMF). We note that the dispersion in the pre- 
dicted properties of the early-type population in the lo- 
cal Universe doe s not vary signi ficantly with other simi- 
lar IMFs such as lSalpeten (|1955l ) - see Table 1 in Kaviraj 
et al. (2005a). 



Notwithstanding the large set of free parameters in- 
herent to any semi-analytical analysis, we note that the 
fiducial values of the parameters in the model are cali- 
brated either through observational data (e.g. the SF ef- 
ficiency) or through numerical simulations (e.g. the merger 
law). Most importantly, the full suite of calibrated pa- 
rameters that drives the model predictions reproduces 
a fundamental set of spectro-photometric properties of 
the observed galaxy population in the Universe. GALICS 
produces good agreement to the galaxy luminosity func- 
tions observed by the 2dF survey in the B and K bands 
(jCross fc 2dF collaboration] l200lh . The (B - V) colours 
of spiral galaxie s closely follow the observed data of 
iButa et alT(f l994) , both in terms of average values and scat- 
ter. Satisfactory fits to the Faber- Jackson relation for early- 
types (and the Tully-Fisher relation for disks) are obtained 
and predictions show good agreement with the early-type 
Fundamental Plane, with the model predicting the observed 
morphological mix of the local Universe with a reasonable 
degree of accuracy. As we describe in detail in the next sec- 
tion, GALICS has been specifically tested against a wide 
variety of multi- wavelength (UV-optical) early-type galaxy 
properties across a wide range in redshift (0 < z < 1.5), 
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making it the ideal tool to map the early-type progenitor 
set over the redshift range < z < 1 . 

We now turn briefly to the impact of the free parame- 
ters on our subsequent analysis and explore the potential 
uncertainty in our results, given the freedom that we may 
have in the values of the parameters in the model. Recalling 
that the primary aim of the paper is to map the photomet- 
ric properties of the early- type progenitor set, we note that 
the parameters that potentially affect the analysis the most 
are (a) those that are responsible for the morphological mix 
of the Universe (since this determines the composition of 
the progenitor set at any given time) and (b) those that 
strongly affect galaxy luminosities/colours (since the prob- 
abilistic prescriptions designed to include late-type progen- 
itors are given as a function of luminosities and colours). 

As described above the key parameter that drives mor- 
phology is the merger law, in particular the transition 
mass ratio (~ 1 : 3) at which morphological transforma- 
tions takes place (i.e. disks are disrupted and bulges form) . 
However, the merger recipe is constrained reasonably ro- 
bustly through numerical simulations and the galaxy pop- 
ulation predicted by the model is consistent with the mor- 
phological mix of observed Universe, indicating that the 
composition of the progenitor set is predicted with an ac- 
ceptable level of accuracy across our target redshift range. 

In a similar vein, galaxy luminosities/colours are driven 
by a suite of parameters including the SF efficiency, SN 
feedback, IMF and dust recipes. This set of parameters is 
constrained through comparison of the fiducial model to the 
properties of observed galaxies. While we have, in principle, 
freedom to change these parameters, calibrations to multi- 
ple observational facts in the real Universe, such as the ones 
described in the previous section, severely reduce this free- 
dom in practice. In other words, if individual parameters 
are altered arbitrarily then the reproduction of the galaxy 
properties in the model would fail. In this sense, the analy- 
sis presented in the paper is stable, in that they are driven 
by values of the free parameters that are based either on 
observations/hydrodynamical simulations and which repro- 
duce the galaxy properties in the present-day Universe with 
a reasonable degree of accuracy. 

2.4. Reproduction of the early-type galaxy population 

Since our aim is to map the photometric properties of early- 
type progenitors, it is important that the model reproduces 
the multi-wavelength properties of the early-type popula- 
tion and their evolution with redshift. A well-calibrated 
model is clearly needed for the progenitor predictions to be 
reliable. Note that in the redshift interval (z < 1) and for 
the wavelength range {B and V band) studied in this pa- 
per, the population synthesis model s employed by the fidu- 
cial GALICS model - STARDUST dDevriendt et alJll999h 
- provide virtuall y identical results to other common ly use d 
models lfYl l2003tlBruzual fc Charioll2003tlMlTrastonll2005h . 

In addition to the reproduction of the general galaxy 
population, the fiducial GALICS model has been specif- 
ically tested in the context of early- type galaxies, across 
virtually the entire redshift range over which early-types 
have been observed (0 < z < 1.5). GALICS accurately re- 
produces the optical colour magnitude relations (CMRs) 
of the early-type population (in dens e environments) an d 
their evolution from z = to z ~ 1.23 (|Kavirai et al.ll2005[) . 
The predicted evolution of both the gradient and the scat- 



ter in the optical CMRs are consistent within errors with 
various observational studies which use a variety of opti- 
cal colours (see Figure 8 in iKavirai et al.l (|2005l )). More 
massive ellipticals are predicted to be older (although they 
assemble more recently) and more metal-ri ch than their 
less m assive counterparts (see Figure 3 in IKavirai et al.l 
(2005))). The star formation histories of elliptical galaxies 
are shown to be quasi-monolithic, which enables the model 
to maintain the correct gradient and scatter over the en- 
tire range in redshift (0 < z < 1.27) at which observa- 
tiona l early-type st u dies h ave been conducted (see Figure 
5 in Kavirai et al.l ((2005)). Note that the decoupling of 
the assembly history of elliptical galaxies from their star 
formation history described abov e has also been found in 



other semi-analytical work (e.g. 


Kauffmann et al.l Il996bt 


Bausch et alJll996t iDe Lucia et al. 


2006). although thev do 



not study the evolution of the optical colours with redshift. 
Finally, GALICS is the only semi-analytic model to have 
been tested against the new generation of UV photomet- 
ric data, shortward of 3000A, made available by the recent 
GALEX mission. Kaviraj et al. (2007) find excellent quanti- 
tative agreement between the predictions of this model and 
the observed t/U-optical photometry of ~ 2100 early-type 
galaxies in the nearby (0 < z < 0.11) Universe that have 
been imaged by both GALEX and the SDSS DR3. 

Given its good reproduction of early-type photometry 
across a wide wavelength range and its successful prediction 
of early-type colour evolution to high redshifts (z ~ 1-23), 
GALICS is a useful and well-calibrated tool with which 
to follow the progenitor set of early-type galaxies and the 
evolution of that progenitor set in the redshift range < 
z < 1. We note that the emphasis in this study is not to 
focus on the properties (e.g. star formation and assembly 
histories) of early-type remnants at z — (which is the 
major thrust of previous semi-analytical studies of elliptical 
galaxies), but to focus on the properties of the progenitor 
galaxies over the redshift range covered by the bulk of the 
recent and forthcoming observational surveys (0 < z < 1). 

3. Dissecting the progenitor set: morphologies of 
progenitors and epochs of last mergers 

Early-type galaxies have an assortment of star formation 
histories (SFHs). Central to this study is the epoch at 
which the last merger, that finally creates the early-type 
remnant, takes place. After this event, the early-type rem- 
nant evolves to present day without further interactions 
with other galaxies. In the discussion that follows, we refer 
to the look-back time to this 'last-merger' as the dynamical 
age of a galaxy. Note that the last mergers typically have 
mass ratios of 1:5 or higher. 

As our subsequent analysis of the progenitor popula- 
tion involves environments of model galaxies, an explana- 
tion of the definition of this quantity is necessary. Galaxy 
environments in the model are driven by the mass of the 
dark matter (DM) halo in which they are embedded. At 
2 = 0, DM halo masses greater than ~ 10 14 M Q corre- 
spond to 'cluster' environments, while halo masses between 
~ 10 13 M Q and ~ 1O 14 M correspond to 'groups'. All other 
halo masses correspond to the 'field'. At higher redshifts 
these definitions do not strictly hold since the DM halo pop- 
ulation is evolving - for example, the largest haloes at z = 1 
are likely to be roug hly half their size at present day (e.g 
Ivan den Boschll2002l ). We take this mass accretion history 
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Table 1. Free parameters that affect the evolution of the baryonic Universe, see Section 2.3 for more details. 



Parameter 


Description 


Fiducial value 


Source of constraints 


P 


Inverse of star formation efficiency 


50 


Kennicutt (1998) 


e 


Inverse of mass loading for feedback 


0.1 


Martin et al. (2002) 


X 


Galaxy merger power law 


3.333 


Numerical simulations e.g. Walker et al. (1996) 



Qb Baryon fraction 0.02 h~ 

S-S merging normalisation 0.017 

? Recycling efficiency 0.3 

k Ratio of burst-to- bulge radius 0.1 
Dust ISM Extinction and emission due to dust 
IMF Initial mass function 



H abundance in QSO absorption lines (T^tler^e^ak, 1996) 
Makino fc Hut (1997) 
Set to fiducial value in Hatton et al. (2003) 
Set to fiducial value in Hatton et al. (2003) 
Guiderdoni et al. (1987), Desert et al. (1990) 
Kennicutt (1983) 



into account when specifying the environments of galaxies 
at high redshif t . The mass accretion history is taken from 
van den Bosch! (|2002l see their Figure 5). 

The number of early-types galaxies in the simulation, 
based on the definitions presented above, is 5418. 279 early- 
types are in clusters (~ 5%), 1081 (~ 20%) are in groups 
and 4058 (~ 75%) are in the field. At z = there are 8 
clusters identified in the GALICS simulation box. 

Figure [2] indicates the last merger redshifts of the sam- 
ple of early-types in the model, split by the environment of 
the remnant at z = 0. Recall that after the last merger, the 
early-type remnant evolves to present day without further 
interactions with other galaxies. The inset in FigureOshows 
histograms of the last merger redshifts, again split by en- 
vironment. As expected, we find that, in all environments, 
larger early-types a re assembled later ( although their stars 
are generally older (iKaviraj et al.l 20051). This re sult is con- 



sistent with the findings of De Lucia et al.l (|2006J) , although 
we note that they did not consider early-types in different 
environments separately. 

Cluster galaxies (at least those brighter than L*) have 
significantly larger dynamical ages - morphological trans- 
formations in clusters therefore proceed more quickly than 
in all other environments. This point is made more clearly 
in Figure [21 where we plot the cumulative fraction of early- 
type galaxies which have already had their last merger. We 
find that, on average, without reference to environment, 
only 35 percent of early- type galaxies are 'in place' (i.e. they 
evolve without further interactions with other galaxies) by 
z = 1 (black line). The rest are still 'in pieces'. In terms 
of morphological transformations, cluster environments are 
special, in that early-type morphologies are attained signifi- 
cantly faster in clusters (red curve) , with almost 70 percent 
of early-type galaxies having undergone their last merger 
by z = 1. 

Before the last merger occurs, the morphology of the 
progenitors is not necessarily early-type. Figure S] shows 
the morphologies of progenitors in binarys mergers as a 
function of redshift. Non-binary mergers do happen but 
are rare, and only take place, in the model, at red- 
shifts greate^_thaii_l^^n agreement with previous work 
(e.g. iKhochfar fc Burkertl 12001 we find that, in the local 
Universe, mergers between early-type progenitors make up 
less than 20 percent of the merger activity. All other merg- 
ers contain at least one spiral progenitor. Mergers involving 
solely spiral progenitors increasingly dominate at higher 
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Fig. 2. Last merger redshifts of the early- types in the 
model, split by environment of the remnant at z = 0. 
INSET: Histograms of last merger redshifts shown in the 
top panel, split by environment. 



redshift and dominate the merger activity beyond z = 1 
(see also lKang et al.ll2007t) . 

Having provided a picture of the merger activity within 
the progenitor set, it is instructive to look at the fraction 
of the progenitor set which is made up of a certain mor- 
phological type as a function of redshift. Figures [5] and |B] 
show the number and mass fractions respectively of progen- 
itors of different morphological types, in the redshift range 
< z < 3 and split by environment and luminosity of the 
early- type remnant at z = 0. We find that, averaging across 
all environments, at z ~ 1, less than 50 percent of the stel- 
lar mass which ends up in early-types today is actually in 
early- type progenitors at this redshift. Faster morphologi- 
cal transformations in cluster environments means that this 
value is ~ 65 percent in clusters at z ~ 1. As a result, look- 
ing only at early-type galaxies at z ~ 1 does not take into 
account almost half the stellar mass in the progenitor set. 
In other words, the mass in the progenitor set doubles be- 
tween z = 1 and z = 0. A similar observational result was 
found by Bell at al. (2004) and Faber et al. (2007), who used 
the optical 'red sequence' as a proxy for the progenitor set 
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Fig. 3. Cumulative fraction of early-type galaxies which 
have had their last merger as a function of redshift. 
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Fig. 5. Number (top) and mass (bottom) fractions con- 
tained in progenitors of different morphological types in 
the redshift range < z < 3, split by the environment of 
the early-type remnant. 



Fig. 4. Morphologies of progenitors in binary mergers as a 
function of redshift. Non-binary mergers do happen but are 
rare, and only take place, in the model, at redshifts greater 
than z ~ 1.5. Mergers at intermediate and high redshift are 
dominated by pairs of progenitors which contain at least 
one spiral progenitor. 



of present-day early-type galaxies (see also the discussion 
in Section 6 below). 

The bias does not arise simply because some progen- 
itor mass is not taken into account, but because the age 
profile of the mass in progenitors of different morphologi- 
cal types tends to vary. We illustrate this point in Figure 
El The top panel shows the average iVt/F-weighted ages 
of progenitors of different morphological types. The NUV 
weighting, generated using the GALEX (Martin el al. 2005) 
NUV filter, is heavily dominated by stars formed within the 
last 0.5 Gyrs of look-back time. At all redshifts, early-type 
progenitors have higher 7V[/l/-wcightcd ages, because the 
mass fraction contributed by recent star formation (RSF) 
i.e. within the last 0.5 Gyrs is smaller than for spiral pro- 
genitors. The differences between elliptical and spiral pro- 
genitors are most pronounced at low redshift. The bottom 
panel shows the fraction of the RSF across the progenitor 



set that is contained in each morphological type. This plot 
has to be interpreted in conjunction with the mass fractions 
hosted by each morphological type as a function of redshift. 
For example, at z ~ 0.1, although spiral progenitors host 
~ 40 percent of the total RSF in the progenitor set, they 
only constitute ~ 30 percent of the mass in the progenitor 
set (see bottom panel of Figure [5]). Early- type progenitors 
(elliptical and SO taken together) contribute ~ 60 percent 
of the RSF - but they also constitute ~ 70 percent of the 
total mass in the progenitor set. Therefore at z ~ 0.1 spi- 
ral progenitors host 1.5 times the amount of RSF per unit 
masss than their early-type counterparts. At higher redshift 
the balance of RSF contained in each morphological type 
moves towards spiral progenitors, partly because they are 
more spirals in the Universe than early- types. 

Figure [7] illustrates that an increasingly larger frac- 
tion of RSF in the progenitor set is contained in late-type 
systems at increasing redshift. In the context of colour- 
magnitude relations (CMRs), which are often used to age- 
date early-type populations at all redshifts, the exclusion 
of spiral progenitors at high redshift biases the CMR to- 
wards redder colours and does not give a proper indication 
of the age of all the stellar mass that eventually constitutes 
present-day early-type galaxies. This is particularly true if 
blue filters (e.g. [/-band or shorter wavelengths) are used 
in the age-dating analysis. 
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Fig. 6. Number (top) and mass (bottom) fractions con- 
tained in progenitors of different morphological types in 
the redshift range < z < 3, split by the luminosity of the 
early-type remnant at present-day. 



4. The spiral progenitors 

One of the aims of this study is to provide a means of in- 
cluding spiral galaxies observed at high redshift which may 
be progenitors into studies of early-type galaxy evolution, 
and thus correct, at least partially, for progenitor bias. We 
therefore focus on spiral progenitors predicted by the model 
and compare their photometric properties to the general 
spiral population. 

4.1. The luminosity function of spiral progenitors 

Wc begin by studying the luminosity function (LF) of spiral 
progenitors. We are interested in studying how the lumi- 
nosities of spiral progenitors compare to the general spiral 
population and what fraction of spirals, at a given lumi- 
nosity, are early-type progenitors. In Figure [8] wc show the 
evolution of the B-band LF of spiral progenitors - we also 
show the spiral progenitors separated by the environment of 
their corresponding early-type remnant at present-day. The 
left hand column illustrates the evolution of the spiral LFs 
- the yellow curve denotes the LF of the general spiral pop- 
ulation (i.e. progenitors + non-progenitors) and the black 
curve the LF of spiral progenitors only. The LFs of spiral 
progenitors whose early- type remnants are, at z = 0, in 
clusters, groups and the field, are shown in red, green and 
blue respectively. The right-hand column shows the frac- 
tion of spiral galaxies which are early-type progenitors as a 
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Fig. 7. TOP: Average AC/T^-weighted ages of progenitors 
of different morphological types. The NUV weighting is 
heavily dominated by stars formed in these progenitors 
within the last 0.5 Gyrs of look-back time. The NUV 
weighting was generated using the GALEX (Martin el al. 
2005) NUV filter. BOTTOM: The mass fraction of recently 
formed stars (age < 1 Gyrs old) across the progenitor set 
which is contained in progenitors of each morphological 
type. Note that the bottom panel has to be interpreted 
in conjunction with the mass fractions hosted by each mor- 
phological type as a function of redshift. For example, at 
z ~ 0.1, although spiral progenitors host ~ 40 percent of 
the total RSF in the progenitor set, they only constitute ~ 
30 percent of the mass in the progenitor set (see bottom 
panel of Figure [3]). Early- type progenitors (elliptical and 
SO taken together) contribute ~ 60 percent of the RSF - 
but they also constitute ~ 70 percent of the total mass in 
the progenitor set. Therefore at z ~ 0.1 spiral progenitors 
host 1.5 times the amount of RSF per unit mass than their 
early- type counterparts. 



function of their B-band luminosities. The colour coding is 
identical to that used for the left-hand column. 

It is apparent that there is a greater preponderance of 
progenitors amongst larger spirals at all redshifts. For ex- 
ample, at low redshifts (z < 0.1), 20 to 40 percent of spirals 
with M(B) < —20.5 are early-type progenitors. At inter- 
mediate redshifts (0.3 < z < 0.52), these values rise to 30 
and 60 percent respectively. At high redshift (z ~ 1) spirals 
with M(B) < —21.5 have more than a 60 percent proba- 
bility of being an early-type progenitor, while spirals with 
-20 < M(B) < -21.5 have between a 30 and 40 percent 
chance of being early-type progenitors. The falling progen- 
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Fig. 8. Evolution of the B-band LF of spiral progenitors. The left hand column illustrates the evolution of the spiral 
LFs - the yellow curve denotes the LF of the general spiral population (i.e. progenitors + non-progenitors) and the black 
curve the LF of spiral progenitors only. The LFs of spiral progenitors whose early-type remnants are, at z = 0, in clusters, 
groups and the field, are shown in red, green and blue respectively. The right-hand column shows the fraction of spiral 
galaxies which are early-type progenitors as a function of their B-band luminosities. The colour coding is identical to 
that used for the left-hand column. 



itor fractions towards lower redshift are partly due to the 
changing morphological mix of the Universe. 

4.2. The colour magnitude space of spiral progenitors 

While investigating the LFs of spiral progenitors is useful in 
indicating the probability that a spiral of a given luminosity 
has of being a progenitor, it is also desirable to explore 
the colour-magnitude (CM) space of the spiral population, 
so that we can separate progenitor spirals better from the 
general population at a given luminosity. 

In Figure Owe compare the (B — V) colours of the gen- 
eral spiral population to the (B — V) colours of spiral pro- 



genitors. The left-hand column shows the spiral (B — V) 
CMR from z = to z = 1. Black dots represent the spi- 
ral galaxies and red dots represent spiral progenitors. In the 
right hand column we show the fraction of spiral progenitors 
across the (B — V) CMR. The fractions are indicated us- 
ing the colour coding shown in the legend. Warmer colours 
indicate a higher progenitor fraction (red implies a progen- 
itor fraction of 1, black represents a progenitor fraction of 
and parts of the CM space without any galaxies are not 
colour-coded). 

At z ~ 0.1, spirals with -21.5 < M(B) < -20.5 have 
~ 30 percent chance of being a progenitor. For larger spi- 
rals, those with red (B — V) colours (i.e. [B — V) > 0.8) 
have ~ 60 percent chance of being a progenitor, while the 
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Fig. 9. (B — V) colours of the general population compared to the (B — V) colours 
of those spirals which are progenitors of present-day early-type galaxies. Black dots 
represent spiral galaxies and red dots represent spiral progenitors. The left-hand 
column shows the spiral {B — V) CMR from z = to z = 1. In the right hand 
column we show the corresponding (binned) fraction of spiral progenitors across 
the (B — V) CM space. The fraction values are indicated using the colour coding 
shown in the legend. Warmer colours indicate a higher progenitor fraction (red 
implies a progenitor fraction of 1, black represents a progenitor fraction of and 
parts of the CM space without any galaxies are not colour-coded). 



corresponding probability for bluer spirals is 30 to 50 per- 
cent. 

At intermediate rcdshift (z ~ 0.5), black spirals, with 
-21.5 < M(B) < -20.5 and (B — V) > 0.6, have - 30 per- 
cent probability of being an early-type progenitor, while 
blue spirals in the same luminosity range have a low pro- 
genitor probability. For larger spirals at these redshifts, 
the probabilities are appreciably higher - red spirals with 
(B — V) > 0.7 have between a 75 and 95 percent chance of 
being progenitors, while 50 to 75 percent of blue spirals in 



this luminosity range are progenitors. The situation at high 
redshift z ~ 1 is similar to that at intermediate rcdshift. 



5. Progenitor evolution in clusters 

Before the advent of large scale surveys, dense regions of the 
Universe were often targetted for e arly-type galaxy studies , 
both at low and high redshift (e.g. iBower et all 



Stanford et all 119981: Ivan Dokkum et all 1199 



2001J) . While studies of dense regions are attractive for a va- 
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riety of reasons (e.g. lEllisI 120021: Ivan Dokkum1l2004h . a key 
benefit is statistical convenience - clusters provide access 
to large homogeneous samples of luminous objects at all 
redshifts. It has been usual to 'connect' results from clus- 
ter studies over large redshift ranges to determine (at least 
qualitatively) the chronology of galaxy evolution. 

In this section, we investigate progenitors of present- 
day cluster early-types, which are themselves in dense re- 
gions at z > 0. The motivation for this investigation is two 
fold. Firstly (and most importantly), it provides a com- 
parison to the vast literature of 'cluster' early- type studies. 
Secondly, this version of the GALICS model has been accu- 
rately calibrated to match the optical CMRs of early-types 
in de nse regions from low to high redshift (jKavirai et al.l 
120051) , which implies that the colours of early- type progen- 
itors (regardless of morphology) are robustly reproduced. 

It is important here to clarify the definitions of 'density' 
that we use to define both 'clusters' at present-day, and 
'dense' regions at high redshift. As mentioned before, model 
'density' is assumed to be a direct function of the mass 
of the DM halo in which galaxies are embedded. At z = 
0, a DM halo mass of 10 14 M Q represents the lower limit 
for a cluster-hosting halo. Observational studies of dense 
regions at high redshift are likely to contain an assortment 
of cluster-type haloes of varying occupancies. Furthermore, 
DM haloes themselves are evolving - on average, the largest 
haloes at z = 1 are likely to be roug hly half their size at 
present day fe.g Ivan den Boschll2002D . To take these two 
points into account, we use a variable lower mass limit for 
'cluster-hosting' haloes at z > 0. At a given redshift this 
lower li mit is calculated from the average mass accretion 
history ([van den Bosch|[20fjl see their Figure 5) applied to 
a 10 14 M o halo. 

Figure [TO] shows the (B - V) CMR in clusters in the 
redshift range < z < 1. Large diamonds indicate pro- 
genitor galaxies - black indicates ellipticals, green indicates 
SOs and blue indicates spiral galaxies. Small crosses indicate 
galaxies which do not contribute to the mass in present-day 
cluster early-types. All early-type galaxies in dense regions 
are, not unexpectedly, progenitors of cluster early-types at 
present day. The top panel in Figure QT] shows the fraction 
of spiral galaxies in dense regions at high redshift (split by 
luminosity) which are progenitors of early- types at z = 0. 
The bottom panel shows the offset in (B — V), with re- 
spect to elliptical progenitors, of the SO and spiral progen- 
itor galaxies. The offsets are shown split by luminosity. 

We find that at high redshift (z ~ 1), up to 40 per- 
cent of large spirals (—23 < M(V) < —21) are progenitors, 
whereas only ~10 percent of small spirals (M(V) > —21) 
are members of the progenitor set. Large spirals are four 
times more likely to be progenitors than small spirals, re- 
gardless of redshift, in the redshift range < z < 1. 
Elliptical galaxies form the reddest locus in (B — V). SO 
galaxies show an average offset of -0.04 compared to the 
elliptical population, regardless of luminosity. Large spiral 
progenitors show an average (B — V) offset of -0.05 com- 
pared to the elliptical population at high redshift, mainly 
because the scatter in the elliptical colours also tends to be 
large at high redshift. At low redshift the offset is more pro- 
nounced - large spiral progenitors are upto 0.1 mags bluer 
in (B — V) than the elliptical population. 
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Fig. 10. The (B - V) CMR in clusters in the redshift range 
< z < 1. Large diamonds indicate the progenitors of 
present-day cluster early-type galaxies. Small crosses in- 
dicate galaxies which do not contribute to the mass in 
present-day cluster early-types. Note that all elliptical and 
SO galaxies in dense regions are, not unexpectedly, progen- 
itors of present-day cluster early-types. 



6. The 'red sequence' as a proxy for the progenitor 

set 

Early-type galaxies in clusters tend to preferentially popu- 
late the reddest parts of the CM space. In this section we 
investigate whether the sample of galaxies (without refer- 
ence to morphology), within the 'red sequence' can be used 
as a proxy for the progenitor set. We define the red sequence 
as the galaxy population which occupies the part of the CM 
space which is dominated by early- type galaxies. In Figures 
[T2]andrr3]thc CM space dominated by early-type galaxies is 
shown in grey - this region is determined by a progressive 
one-sigma fit to the colours of the early-type population. 
It therefore contains, on average, 68 percent of the early- 
type population within it. Large diamonds indicate galax- 
ies which are part of the progenitor set. Galaxies which 
are not part of the progenitor set are shown using small 
dots. Galaxies in the red sequence are circled. It is appar- 
ent that the red sequence misses blue progenitor galaxies, 
both early-type and late-type. Figure [T2] shows the evolu- 
tion of the red sequence population from low to interme- 
diate redshift and Figure [13] the corresponding evolution 
from intermediate to high redshift. 

In Figure Q3] we compare galaxies in the actual progen- 
itor set to those in the red sequence. The top panel shows 
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Fig. 11. TOP PANEL: The fraction of spiral galaxies in 
dense regions (split by luminosity) which are progenitors of 
cluster early-types at z = 0. Note that all elliptical and SO 
galaxies in dense regions are, not unexpectedly, progenitors 
of present-day cluster early-types. BOTTOM PANEL: The 
offset in average (B — V) colour, with respect to elliptical 
progenitors, of SO and spiral progenitors in dense regions. 
The offsets are shown split by luminosity 



the number ratio between the progenitor set population and 
the red sequence population, split by rcdshift and luminos- 
ity, while the bottom panel shows the mass ratio between 
the progenitor set population and the red sequence popu- 
lation split by redshift and luminosity. It is apparent that 
large galaxies (—23 < M(V) < —21) in the red sequence 
trace the progenitor set well in terms of number and mass 
but that the relationship breaks down as wc move towards 
the lower end of the luminosity function (M(V) > —21). 
The number and mass fractions remain stable, as a func- 
tion of the luminosity slices shown in Figure [T4l within the 
redshift range explored in this study (0 < z < 1). 

Luminosity evolution studies which use the red sequence 
as a proxy for the early-type population can therefore 
achieve accurate results only for the massive end of the lu- 
minosity function. Note that the results of two such recent 
studies, Bell et al. (2004) and Faber et al. (2007), are accu- 
rate since their analysis of the luminosity evolution of red 
galaxies is indeed restricted to massive galaxies on the red 
sequence (e.g. M(V) < -21 in Bell et al. 2004). However, 
it should be noted that the red sequence is less reliable 
as a proxy for the progenitor set further down the lumi- 
nosity function. In addition, the red sequence, almost by 
definition, misses contributions due to early-types which 



Fig. 12. The composition of the 'red sequence', defined as 
the galaxy population which occupies the part of the CM 
space dominated by early-type galaxies (shown in grey), 
compared with the progenitor set in cluster populations. 
The diamonds indicate galaxies which are part of the pro- 
genitor set. Galaxies which are not part of the progenitor 
set are shown using small crosses. Galaxies in the red se- 
quence are circled. The 'red sequence' misses blue galaxies, 
both early-type and late-type. The top panel shows the 'red 
sequence' at z = 0.11, while the bottom panel shows the 
'red sequence' at z = 0.32. 



lie blueward of it. This is increasingly true at higher red- 
shift - hence conclusions based on the colours of the red 
sequence should be applied carefully to early-type evolu- 
tion as deeper or higher redshift data becomes available in 
the future. 

7. Summary and discussion 

We have presented a comprehensive theoretical study, us- 
ing a realistic semi-analytical model in the framework of 
the standard LCDM paradigm, of the photometric proper- 
ties of the progenitors of present-day early- type galaxies. 
We have mapped, in detail, the properties of these progen- 
itors, regardless of their morphology, in the redshift range 
< z < 1 (where the bulk of current large-scale surveys are 
focussed), as a function of the luminosity and environment 
of the early-type remnant at present day. We have also de- 
veloped probabilistic prescriptions which provide a means 
of including spiral (i.e. non early-type) progenitors at in- 
termediate and high redshifts, based on their luminosity 
and optical (BVK) colours. The prescriptions developed 
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Fig. 13. The composition of the 'red sequence', defined as 
the galaxy population which occupies the part of the CM 
space dominated by early- type galaxies (shown in grey), 
compared with the progenitor set in cluster populations. 
The diamonds indicate galaxies which are part of the pro- 
genitor set. Galaxies which are not part of the progenitor set 
are shown using small crosses. Galaxies in the red sequence 
are circled. The 'red sequence' misses blue galaxies, both 
early-type and late-type. This plot shows the evolution of 
the 'red sequence' from intermediate to high redshift. The 
top panel shows the 'red sequence' at z = 0.57, while the 
bottom panel shows the 'red sequence' at z = 1. 



here can potentially be used to address, from the perspec- 
tive of the standard modelss, the issue of 'progenitor bias', 
whereby the exclusion of late-type progenitors in observa- 
tional studies can lead to inaccurate conclusions regarding 
the evolution of the early-type population over cosmic time. 
While simple prescriptions designed to alleviate progenitor 
bias do exist and have been used in some observational 
studies, they do not reflect the mass assembly and morpho- 
logical transformations inherent to the standard model. 

Against the backdrop of the impending or recent release 
of data from large-scale surveys at high redshift, the results 
of this study are hopefully timely, providing a picture of 
how the mass in present-day early-types, across a range of 
luminosities and environments, have assembled and explor- 
ing the extent to which spiral progenitors play a role in 
this process. Our main conclusions can be summarised as 
follows: 

— Larger early-types in all environments arc assem- 
bled later than their less massive counterparts. 



Fig. 14. Comparison between galaxies in the actual progen- 
itor set to those in the red sequence. TOP: Number ratio 
between the progenitor set population and the red sequence 
population split by redshift and luminosity. BOTTOM: 
Mass ratio between the progenitor set population and the 
red sequence population split by redshift and luminosity. It 
is apparent that large galaxies in the red sequence trace the 
progenitor set well in terms of number and mass but that 
the relationship breaks down rather rapidly as we move to- 
wards the lower end of the luminosity function. 



However, their stellar populations are generally older. 
This result is consistent with previous studies of 
early- type galaxies in t he semi -a nalytical framework 



(e.g. Kauffmann et al. 1996bt Baueh et al.l 119961 : 



iKavirai et al.ll2005HDe Lucia et al. 2006 ) 



On average, without reference to environment, only 35 
percent of early-type galaxies are in place by z = 1 
and evolve without further interactions with other 
galaxies thereafter. Morphological transformations are 
significantly faster in cluster environments (where the 
vast majority of early- type studies have been based 
before the advent of large-scale surveys). In clusters 
almost 70 percent of early- types are in place by z = 1. 
In other words, the probability of a 'major merger', 
which creates an early-type remnant is low after z = 1 
in cluster type environments. 

Averaging across all environments, at z ~ 1, less than 
50 percent of the stellar mass which ends up in early- 
types today is actually in early-type progenitors at this 
redshift. This value is around 65 percent in clusters 
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owing to faster morphological transformation in this 
environment. In other words, looking only at early-type 
progenitors does not take into account almost half the 
mass in the progenitor set - the progenitor set doubles 
in mass in the redshift range < z < 1. 

— Progenitor bias does not arise simply because (late- 
type) progenitor mass is missed, but also because 
the age profile of mass in progenitors of different 
morphological types tend to vary. Spiral progenitors 
are typically 'bluer' (at a given redshift) because 
they host more recently formed stars than early-type 
progenitors. Hence, age-dating the progenitor set using 
an 'early-type only' CMR, i.e. after excluding the spiral 
progenitors, biases the CMR towards redder colours 
and overestimates the average age of the progenitor set. 

— One of the aims of this study is to provide a means of 
including spiral progenitors into studies of early-type 
evolution, and thus correct, at least partially, for 
progenitor bias. We have therefore focussed on spiral 
progenitors in the model and compared their properties, 
in detail, to the general spiral population. 

— There is a greater preponderance of progenitors 
amongst larger spirals at all rcdshifts. At low red- 
shifts (z < 0.1), 20 to 40 percent of spirals with 
M(B) < —20.5 are early-type progenitors. At interme- 
diate redshifts 0.3 < z < 0.5, these values rise to 30 and 
60 percent respectively. At high redshift (z ~ 1) spirals 
with M(B) < —21.5 have more than a 60 percent 
probability of being an early-type progenitor while 
spirals with -20 < M(B) < -21.5 have between a 30 
and 40 percent chance of being early-type progenitors. 
The falling progenitor fractions towards lower redshift 
are partly due to the changing morphological mix of 
the Universe. 

— The colour-magnitude space of the spiral popu- 
lation provides a slightly better route to identi- 
fying spiral progenitors. At z ~ 0.1, spirals with 
-21.5 < M(B) < -20.5 have ~ 30 percent chance 
of being a progenitor. For larger spirals, those with 
red (B - V) colours, i.e. (B - V) > 0.8, have - 60 
percent chance of being a progenitor, while the corre- 
sponding probability for bluer spirals is 30 to 50 percent. 

At intermediate redshift (z ~ 0.5), black spirals, with 
-21.5 < M(B) < -20.5 and (B — V) > 0.6, have - 30 
percent probability of being an early-type progenitor, 
while blue spirals in the same luminosity range have a 
low progenitor probability. For larger spirals at these 
redshifts, the probabilities are appreciably higher - red 
spirals with (B — V) > 0.7 have between a 75 and 
95 percent chance of being progenitors, while 50 to 
75 percent of blue spirals in this luminosity range are 
progenitors. The situation at high redshift z ~ 1 is 
similar to that at intermediate redshift. The trends in 
the (V — K) colour are similar to those in (B — V). 

— We suggest one example of how the probabilites derived 
in Section 4 could be used to 'include' spiral progenitors 
at high redshifts (z > 0.5). Modern optical surveys 
(e.g. MUSYC) can be used to trace the rest-frame 



UV photometry of galaxies at high redshifts. Given 
its sensitivity to young stars, the rest-frame UV is an 
excellent photometric tracer of recent star formation, 
i.e. the mass fraction of stars formed within the last 
Gyr. Clearly, if one had access to the rest-frame UV 
over a range of redshifts, one could progressively trace 
the build-up of stellar mass in a particular class of 
galaxy over time. Thus, the stellar mass build-up in the 
early-type progenitor set can be approximated by the 
sum of the recent star formation in early-type galaxies 
within a given redshift range plus the recent star 
formation in spiral galaxies weighted by their relevant 
probabilities of being progenitors derived in Section 4. 

— We note that the analysis presented in this paper 
strongly implies that model and observed early-type 
populations should onlys be compared once they have 
been split by environment, since the assembly history 
of early- type galaxies is a function of environment. 
For example, the combined model population (at 
any redshift) should not be compared to an observed 
population which may be drawn predominantly from 
rich groups and clusters and vice- versa. 

— Finally we have explored the correspondence between 
the progenitor set and the 'red sequence', defined as the 
part of the CM parameter space which is dominated 
by early-type galaxies. We find that galaxies, both 
late and early-type, that fall in this parameter space 
do not necessarily trace the progenitor set well. Large 
galaxies (—23 < M(V) < —21) in the red sequence 
correspond to the progenitor set reasonably well in 
terms of number and mass but the relationship breaks 
down as we go towards the lower end of the luminosity 
function (M(V) > —21). Hence, luminosity evolution 
studies which use the red sequence as a proxy for 
the early-type population, therefore achieve accurate 
results only for the upper end of the luminosity 
function. In addition, the red sequence, almost by defi- 
nition, misses contributions due to early-types which 
lie blueward of it - hence conclusions based on the 
colours of the red sequence should not generally be ap- 
plied to early-type evolution, especially at high redshift. 
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